GMRES/CR and Arnoldi/Lanczos as Matrix Approximation Problems
نویسندگان
چکیده
The GMRES and Arnoldi algorithms, which reduce to the CR and Lanczos algorithms in the symmetric case, both minimize p(A)b over polynomials p of degree n. The difference is that p is normalized at z 0 for GMRES and at z x for Arnoldi. Analogous "ideal GMRES" and "ideal Arnoldi" problems are obtained if one removes b from the discussion and minimizes p(/l)II instead. Investigation of these true and ideal approximation problems gives insight into how fast GMRES converges and how the Arnoldi iteration locates eigenvalues.
منابع مشابه
Tikhonov regularization via flexible Arnoldi reduction
Flexible GMRES, introduced by Saad, is a generalization of the standard GMRES method for the solution of large linear systems of equations. It is based on the flexible Arnoldi process for reducing a large square matrix to a small matrix. We describe how the flexible Arnoldi process can be applied to implement one-parameter and multi-parameter Tikhonov regularization of linear discrete ill-posed...
متن کاملInvestigating the Effects of MINRES with Local Reorthogonalization
The minimum residual method, MINRES, is an iterative method for solving a n x n linear system, Ax = b, where is A is a symmetric matrix. It searches for a vector, xk, in the k th Krylov subspace that minimizes the residual, rk = b − Axk. Due to the symmetric nature, the vectors are computed to be the Lanczos vectors, which simplifies to a three-term recurrence, where vk+1 is computed as a linea...
متن کاملKrylov Subspace Approximation for TEM Simulation in the Time Domain
Forward transient electromagnetic modeling requires the numerical solution of a linear constant-coefficient initial-value problem for the quasi-static Maxwell equations. After discretization in space this problem reduces to a large system of ordinary differential equations, which is typically solved using finite-difference time-stepping. We compare standard time-stepping schemes such as the exp...
متن کاملOn Best Approximations of Polynomials in Matrices in the Matrix 2-Norm
We show that certain matrix approximation problems in the matrix 2-norm have uniquely defined solutions, despite the lack of strict convexity of the matrix 2-norm. The problems we consider are generalizations of the ideal Arnoldi and ideal GMRES approximation problems introduced by Greenbaum and Trefethen [SIAM J. Sci. Comput., 15 (1994), pp. 359–368]. We also discuss general characterizations ...
متن کاملCSSP 13:733-758, 1994 Krylov Space Methods on State-Space Control Models
We give an overview of various Lanczos/Krylov space methods and how they are being used for solving certain problems in Control Systems Theory based on state-space models. The matrix methods used are based on Krylov sequences and are closely related to modern iterative methods for standard matrix problems such as sets of linear equations and eigenvalue calculations. We show how these methods ca...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- SIAM J. Scientific Computing
دوره 15 شماره
صفحات -
تاریخ انتشار 1994